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Abstract: Consider the 2n-by-2n matrix M = imij)f"j^^ 
with mjj = 1 for /, j satisfying |2; — 2n — 1 1 + |2y — 2n — 1 1 < 
2n and m,j — for all other i,j, consisting of a central 
diamond of I's surrounded by O's. When n > 4, the A- 
determinant of the matrix M (as introduced by Robbins and 
Rumsey [7]) is not well-defined. However, if we replace 
the O's by f's, we get a matrix whose A-determinant is well- 
defined and is a polynomial in A and f . The limit of this poly- 
nomial as f ^ is a polynomial in A whose value at A = 1 is 
the number of domino tilings of a 2n-by-2n square. 



1. Lambda-determinants ... 

In section 5 of their article [7], David Robbins and Howard 
Rumsey, Jr. defined a generalization of the determinant of a 
matrix, which they dubbed the A-determinant. It is a ratio- 
nal function of the entries of the matrix along with an extra 
parameter, A ; when A is set equal to — 1 , one obtains the or- 
dinary determinant of the matrix, at least in the case where all 
matrix entries are non-zero. 

In this article 1 will consider certain matrices with many 
vanishing entries. For these matrices, one cannot apply Rob- 
bins and Rumsey's definition literally, but there is still a nat- 
ural way to attempt to compute the A-determinant, by replac- 
ing the zeroes by an indeterminate t, taking the A-determinant 
of the f -perturbed matrix, and then taking the limit as f ^ 
0. In particular, 1 will give a one-parameter family of 0, 1- 
matrices whose nth member is a 2;i-by-2n matrix whose A- 
determinant, defined by this continuity method and then spe- 
cialized to A = 1, is the number of domino-tilings of a 2n-by- 
2n square. 

Let us start by recalling Robbins and Rumsey's first, recur- 
sive definition of the A-determinant. If M is a 1-by-l matrix, 
its A-determinant is its sole entry. Now suppose M is an n-by- 
n matrix, with n > I. Let M^w, M^e, M^w, and Mse denote 
the A -determinants of the (n — l)-by-(n — 1) connected sub- 
matrices in the northwest, northeast, southwest, and southeast 
corners of M, and let Mc be the A-determinant of the centi'al 
connected (n — 2)-by-(n — 2) submatrix of M (we take Mc = 1 
in the case n — 2). As long as Mc is non-zero, we define the 
A-determinant of M as 



det;tM = (MnwMse + IMneMsw)/Mc. 



(1) 



Using this definition, we can calculate the A-determinant of 
the matrix 



1 1 
1 1 




as((l+A)(l+A)+A(l+A)(l+A))/(l) = (l+A)3. Aswas 
pointed out by Robbins and Rumsey (and is easy to check by 
induction), the A-determinant of the n-by-n all-ones matrix is 
(l+A)«("+i)/2. 

When one is using (1) to calculate A -determinants, a sub- 
tle distinction becomes important, namely, the distinction be- 
tween working in Q(A) throughout the recursion and substi- 
tuting a particular value of A at the end (on the one hand), 
and using that particular value of A when performing the re- 
cursion (on the other). Consider, for instance, the 4-by-4 ma- 
trix whose entries are all I's. Its A-determinant is (1 + A)^, 
and if we put A — —1, we get 0. However, if we were to 
use A = — 1 in carrying out the recurrence, we would run into 
trouble, since Mf^w, Mj^e, Msw^ M^e, and Mq all vanish; the 
(—1) -determinant of the matrix is given by the indeterminate 
expression ((0)(0) — (0)(0))/0. If our goal is to be make sense 
of the A-determinant over as broad a class of matrices as pos- 
sible, clearly we should work in Q(A) whenever we can. 

Robbins and Rumsey give another, non-recursive formula 
for the A-determinant of an n-by-n matrix M = {mij}" ■^: 



det^iM)= ^ A^(^)(1+A)^(^'m^ 



(2) 



B £.<:/„ 



Here £/„ is the set of n-by-n alternating-sign matrices B = 
{bi.j)1 -i, P{-) and A^(-) are integer-valued functions on £/„, 
and 



M" 






bi.j 



In more detail: An alternating-sign matrix is a matrix of + 1 's, 
— I's, and O's such that in each row and column, the non- 
zero entries alternate in sign, beginning and ending with a + 1 
(which may be the same entry). If B is an n-by-n alternating- 
sign matrix, we define its inversion number as I{B) ~ Y-bijbr.s 
where the sum is over all 1 < /, j, r,s <n with / < r and j > s 
(if B is a permutation matrix, this coincides with the ordinary 
inversion statistic). We also define N{B) as the number of 
negative entries in B, and P{B) as I{B) —N{B). 

The summation formula (1) for det;^-'^ has exactly the same 
domain of applicability as the recursive formula (2), if we 
work over Q(A) (and if the matrix entries themselves do not 
depend on A); specifically, both formulas apply and give the 
same answer as long as the entries of M that are in the central 
(n — 2)-by-(n — 2) submatrix are all non-zero (these are pre- 
cisely the positions in an n-by-n alternating-sign matrix where 
the entry —1 can occur). However, if we use specific values 
of A in the recursion, the recursive formula can run into prob- 
lems where the summation formula does not. For instance, 
consider the 4-by-4 matrix with all entries equal to 1 ; as we 
have seen, if we try to compute its ( — 1 ) -determinant using the 



specialization of the recurrence to A = — 1, we get an indeter- 
minate result, whereas the formula in terms of alternating-sign 
matrices gives 0. 

Unfortunately, neither of Robbins and Rumsey's two def- 
initions works when central entries of M are equal to zero. 
Nor is a straightforward appeal to continuity going to help us 
to define det;^ M for every M. Consider for instance the family 
of matrices 

(t t t 
t t^/c t 
t t t ^ 

For c and t non-zero, det;i,Mc(f) — (cA + cA^) + (2 A + 
2A^)/^ + (1/c + A^/c)f^, which converges to cX + cX'^ as 
f ^ 0. This limit depends on the value c. Hence, if we at- 
tempt to define the A -determinant of the three-by-three all- 
zeroes matrix by taking a trajectory through that matrix in the 
space of three-by-three matrices and invoking continuity, the 
limit will depend on the trajectory we choose, and may even 
fail to exist. 

Clearly the principled thing to do would be to study conti- 
nuity properties of the A -determinant, and I hope others will 
adopt this approach and undertake a more systematic study of 
what happens when different trajectories through a bad ma- 
trix are taken; this may have some bearing on the issue of how 
Dodgson condensation can be applied to matrices with van- 
ishing connected minors (i.e., with singular connected sub- 
matrices). However, in this article I will take an easier path 
and restrict attention to the trajectory through M in a partic- 
ular direction. Specifically, I will replace the zeroes in M by 
a new variable, t, and see what happens to the A-determinant 
of the perturbed matrix (the "f -perturbation of M") as f ^ 0. 
This may be unprincipled, but it is easy to compute. More- 
over, for many matrices M, the resulting rational function in t 
is not just a Laurent polynomial in t (as it must be, because of 
the summation formula), but is in fact an ordinary polynomial 
in t. In this case, the A-determinant of M can be defined as the 
constant term of this polynomial. 

For example, consider the eight-by-eight matrix 

/O 1 1 0\ 

11110 

1111110 

11111111 

11111111 

1111110 

11110 
VO 1 1 0/ 

It has interior zeroes, so its A-determinant cannot be computed 
by the original formulas (1) and (2). However, if we replace 
all O's by t, we get a polynomial in A and t with 191 terms. 
Replacing t by 0, we get a polynomial in A with a mere 17 
terms, which evaluates to 12988816 when we set A = 1. 

2. ... and Domino-Tilings 

12988816 is also the number of ways to cover an 8-by-8 
square with 32 l-by-2 rectangles, commonly known as domi- 
nos. More generally, the number of ways to cover a 2n-by-2n 



square with 2n^ dominos was computed by Temperley and 
Fisher [9] and (simultaneously) by Kasteleyn [4], and is given 
by the double product 



nn 

j=\k=i 



4cos'^ 



2n + l 



-4cos^ 



%k 
2n + l 



It turns out that if one considers the sub-region of the 2n-by- 
2« square that consists of the 2 central cell in the first and last 
rows, the 4 central cells in the second and second-from-last 
rows, the 6 central cells in the third and third-from-last rows, 
etc., one gets a region whose domino-tilings are enumerated 
by a much simpler expression, namely 

2n(n+l)/2 

This region is called the Aztec diamond of order 8, and was 
first studied in detail in [3], although earlier occurrences of 
the shape appear in the literature. As an example, consider the 
case n = 4. If we associate the 64 cells of the 8-by-8 square 
with the entries of the 8-by-8 matrix considered at the end of 
the previous section, then the entries that contain I's corre- 
spond to the cells that belong to the Aztec diamond of order 
4. 

It turns out that the domino-tilings of the Aztec diamond 
of order n are intimately related to the A-determinant of the 
generic «-by-n matrix. (If the meaning of "generic" is unclear, 
then the reader should imagine that we are working over the 
ring Q{X,a,b,c,...), where a,b,c,... are the entries of the 
matrix; then the entries of the matrix are perforce non-zero, so 
the summation formula for the A-determinant is unproblem- 
atical in this context.) Specifically, if we apply the distribu- 
tive rule to the Robbins-Rumsey summation formula for det;^, 
and ignore the fact that multiplication is commutative, each 
term A^(^'(l + A)'^(^'M^ expands to a sum of 2'^(^' mono- 
mials of the form A'^M^. The resulting sum has 2"^"^^''^ 
terms, and these terms can be put into 1-1 correspondence 
with the 2"("+'V 2 domino-tilings of the Aztec diamond of or- 
der n. (See [3] for details.) 

If it seems paradoxical that the A-determinant of a square 
matrix of I's counts domino-tilings of an Aztec diamond, 
whereas the A-determinant of the 0,1-matrix whose I's form 
an Aztec diamond counts domino-tilings of a square, it may be 
helpful to imagine rotating the Aztec diamond by 45 degrees. 
Here is one of the 64 domino-tilings of the Aztec diamond of 
order 3: 




We can replace the tiling problem by the dual matching prob- 
lem. We define an Aztec diamond graph whose vertices corre- 
spond to the cells of the Aztec diamond, with an edge between 
two vertices when the two corresponding cells are adjacent. 
Then a tiling of the Aztec diamond of order n corresponds to 
a perfect matching of the Aztec diamond graph, that is, a col- 
lection of edges with the property that each vertex of the graph 
belongs to exactly one of the edges in the collection. 

To illustrate, here is the Aztec diamond graph of order 3: 



Aztec diamond graph of order 5 (the edges shown get weight 
1, and the rest get weight 0): 




And here is the perfect matching of the Aztec diamond graph 
of order 3 that corresponds to the domino-tiling of the Aztec 
diamond shown in the first figure: 




However, what was missing from that account was an expla- 
nation of the link with A -determinants. 

A good way to understand this link, without using all the 
machinery of generalized domino-shuffling, is to directly rely 
on Kuo's method of graphical condensation [5]. Let G be a 
weighted Aztec diamond graph of order n. Let Gnw be the 
weighted Aztec diamond graph of order n — 1 derived from G 
by elminating the southernmost n vertices, the southernmost 
2n edges, the easternmost n vertices, and the easternmost 2n 
edges. Define One, Gswi and G^e analogously. Define Gc to 
be the weighted Aztec diamond graph of order n — 2 derived 
from G by eliminating all of the aforementioned vertices and 
edges. Lastly, define wnw to be the weight of the northwest- 
ernmost edge, and define wne, wswi and wse analogously. 
Then Kuo's formula can be stated as 



W(G)W{Gc) == 

WNEWSwW {GnwW {Gse) + WNWWSeW {GtiE)W {Gsw), 



Looking back at the original figure, note that if the northwest, 
northeast, southwest, and southeast tiles are removed, what's 
left is a domino-tiling of the 4-by-4 square. More generally, 
there is a one-to-one correspondence between domino-tilings 
of the 2n-by-2n square and domino-tilings of the Aztec dia- 
mond of order 2« in which there are n{n — l)/2 forced tiles 
in each of the four corners (slanting from southwest to north- 
east in the northwest and southeast corners, and slanting from 
northwest to southeast in the southwest and northeast corners). 
The latter in turn correspond to perfect matchings of the Aztec 
diamond graph of order 2n in which there are n(n — l)/2 
forced edges in each corner. 

One way we can force some edges to be present is to force 
other edges to be absent, and one way to force edges to be ab- 
sent is to work in the setting of weighted enumeration. Given 
an assignment of non-negative weights to the edges of a graph, 
we define the weight of an individual perfect matching as the 
product of the weights of its edges. In the article [6], I showed 
how the method of [3] could be adapted to the general problem 
of finding the sum of the weights of all the perfect matchings 
of an edge-weighted Aztec diamond graph. I applied this to 
the case of 2«-by-2« squares, setting some edge-weights equal 
to 1 and the rest equal to 0, in the fashion shown below for an 



where W(-) denotes the sum of the weights of the perfect 
matchings of the graph in question. In our application, the 
edge-weights wnw, wne, wswi and wse are all O's and I's. 
This formula becomes a recurrence relation if one divides both 
sides by W(Gc) (though this is only sensible if W(Gc) is non- 
zero). If one runs the recurrence for the case at hand, the pref- 
actors wneWsw and wnwWse are always equal to L Hence the 
recurrence takes the simplified form 

W{G) = {W{GnwW{Gse)+W{Gne)W{Gsw))IW{Gc), 

which is equation (1) in the special case X = \. This gives us 
a combinatorial proof of the main claim of this paper, namely, 
that the number of domino-tilings of the 2n-by-2n square can 
be calculated as the (+ 1 ) -determinant of the 2n-by-2n matrix 
whose /, jth entry (for 1 < /, j < 2n) is 1 if |2! — 2m — 1 1 + |2y — 
2n — 1 1 < 2n and is otherwise. 

Kuo's method also gives us a combinatorial interpretation 
to all the numbers that occur in the course of evaluating the 
(+l)-determinantof the matrix (by applying (1) with X = 1): 
these numbers count domino-tilings of regions obtained from 
the Aztec diamond by eliminating the cells that lie in certain 
bands bordering the boundaries of the region. 

For instance, consider the process of computing the (+1)- 



determinant of the 4-by-4 matrix 



M = 




The quantities that turn up in the recursive appHcation of (1) 
are precisely the (+1) -determinants of all the connected sub- 
matrices of M. M itself does double-duty as the matrix of 
(+1) -determinants of all the 1-by-l submatrices of M. The 
(+1) -determinants of the connected 2-by-2 submatrices of M 
form the 3-by-3 matrix 




whose entries are given by the A -condensation formula (1) in 
the special case A = +1. Turning the crank again gives the 
2-by-2 matrix 



any alternating-sign matrix A = («(,;); ';ii ^"d sums those en- 
tries Uij for which |/ — n — 1 1 + |y — « — 1 1 <n, one gets a 
non-negative sum. This allows one to use the f-perturbation 
trick to extend the definition of the A -determinant to the odd- 
by-odd matrix M = (m,y)^"^| with m,- ,■ = 1 for those ij for 
which |/ — n — 1 1 + I J — « — 1 1 <n and niij = for all other 
ij. (In fact, the (+l)-determinant of this matrix is also equal 
to the number of domino-tilings of the 2«-by-2n square. This 
was proved by Trevor Bass and Kezia Charles [1].) 

Moreover, if one takes the intersection of either the odd- 
by-odd or even-by-even diamond pattern with any connected 
square submatrix of M, one gets a smaller pattern which also 
has the property that sum of the corresponding entries of an 
alternating-sign matrix must be non-negative. It would be 
interesting to have a classification of those partial sums of 
the entries of an «-by-« matrix that are non-negative for all 
choices of an alternating-sign matrix. It would also be desir- 
able to extend the results of this paper to the context of the 
octahedron recurrence with general initial conditions, as con- 
sidered in [8]. 



6 6 
6 6 



and turning it a final time gives the 1-by-l matrix 

(36) 

whose sole entry is the number of tilings of the 4-by-4 square. 
The 6's count the domino-tilings of the region obtained from 
the 4-by-4 square by removing three cells from one corner and 
three cells from an adjoining corner 

This combinatorial interpretation of the (+l)-determinants 
of the connected submatrices of M makes it clear that Rob- 
bins and Rumsey's original recursive formula for the A- 
determinant can be applied to all of the matrices in our one- 
parameter family, provided we interpret the indeterminate ex- 
pression 0/0 as whenever it crops up; each such occurrence 
corresponds to a subgraph of the weighted Aztec diamond 
graph whose perfect matchings all have weight 0. 

The one thing missing from this explanation is an explicit 
discussion of the behavior of the (+l)-determinant of the t- 
perturbed matrix. We need to know that it is a polynomial in 
t. But this follows from the fact that it is equal to the sum 
of the weights of all the perfect matchings of the weighted 
graph. Indeed, all the rational functions of t that occur during 
the recursion are polynomials in f, for the same reason. 

This takes care of the case A = +1: the inclusion of t has 
solved the indeterminacy problem. It follows a fortiori that 
for generic A, inclusion of t also lets us carry out the recur- 
rence (1) without encountering indeterminacy. (However, if 
one puts A = — 1, one still encounters indeterminacy, on ac- 
count of the cancellations that occur.) 

One consequence of the main result of this paper is that if 
one takes any 2n-by-2« alternating-sign matrix A ~ {<^i.j)f'i=i 
and sums those entries ajj for which \2i 



2n-l\ + \2j- 
2n — 1 1 < 2n, one gets a non-negative sum. A similar claim 
holds for odd-by-odd alternating-sign matrices: if one takes 
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